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Abstract 

The flowfield of an airfoil oscillated periodically over a 
wide range of reduced frequencies, 0^k^l.6, is studied 
experimentally at chord Reynolds numbers of R c =' 22,000 
and 44,000. The NACA0012 airfoil is pitched sinusoidally 
about one quarter chord between angles of attack (a) of 5° 
and 25°. Detailed flow visualization and phase averaged 
vorticity measurements are carried out for k = 0.2 to docu- 
ment the evolution and the shedding of the dynamic stall 
vortex (DSV). In addition to the DSV, an intense vortex of 
opposite sign originates from the trailing edge just when the 
DSV is shed. After being shed into the wake, the two to- 
gether take the shape of a large M mushroom” while being 
convected away from the airfoil. The unsteady circulation 
around the airfoil and, therefore, the time varying compo- 
nent of the lift is estimated in a novel way from the shed 
vorticity flux and is found to be in good agreement with the 
lift variation reported by others. The delay in the shedding 
of the DSV with increasing k, as observed by previous 
researchers, is documented for the full range of k. The 
DSV, for example, is shed nearly at the maximum a of 25° 
at k = 0.2, but is shed at the minimum a of 5° at k = 
0.8. At low k, the flowfield appears quasi -steady and the 
bluff body shedding corresponding to the maximum a (25°) 
dominates the unsteady fluctuations in the wake. 

1. Introduction 

The phenomenon of dynamic stall on airfoils and lifting 
surfaces in unsteady flow environments has been studied 
experimentally and computationally for many years, both as 
an important practical problem and a challenging fundamen- 
tal one as well. The phenomenon appears on helicopter 
rotor blades, rapidly maneuvering aircraft, fluttering com- 
pressor blades, wind turbines and even fish tails and insect 
wings. It is now well known that the unsteady fluid me- 
chanics of an airfoil pitched sufficiently above the static 
stall limit, is characterized by the formation of a strong 
vortex on the suction surface, known as the dynamic stall 
vortex (abbreviated here as DSV), which is eventually shed 


into the wake. Presence of the DSV on the airfoil upper 
surface causes a dramatic increase in the airfoil lift which, 
however, decreases suddenly when the DSV is shed. Excel- 
lent overviews of this dynamic stall phenomenon can be 
found in the papers by McCroskey 1 and Carr 2 . 

The flowfield of an oscillating airfoil in a uniform 
stream is governed by two time scales: one imposed by the 
pitching motion, T p = l/(2icf); and the other by the free 
stream Velocity and the airfoil chord (conventionally, half 
the chord is used), T fg = c/2U„. Here, f is the oscillation 
frequency, c the airfoil chord and U„ the freestream veloci- 
ty. The reduced frequency, defined as k = 7cfc/U„ thus 
represents the ratio of the two time scales, T fg /T p . For very 
low k, when T p »T fg unsteady effects are relatively weak and 
the flowfield is expected to respond in a quasi -steady man- 
ner. At very high k, on the other hand, intuitively one 
might expect that the outer flow would be unable to res- 
pond to the pitching motion and the moving airfoil might 
appear similar to a bluff body to the flow. 

Most of the previous experimental (e.g. Carr et al. 3 ; 
Leishman 4 ) and computational (e.g., Mehta 5 , Wu et al. 6 ) 
studies on the phenomenon are confined to low values of k 
( < 0.4). Small amplitude pitching motion at very high k has 
been studied by Koochesfahani 7 , among others, but the 
small amplitude (±4°) apparently precluded the possibility 
of dynamic stall. A low Reynolds number experimental 
study by Ohmi et al. 8 covered a higher range of k but the 
effect studied was more complex, viz., that due to a com- 
bined oscillating and translating motion of the airfoil. The 
latter study showed that at relatively low k (= 0.63) the 
flowfield is primarily governed by the translatory motion, 
but at high k (= 3.2 and 6.6) the oscillatory motion con- 
trols the flowfield. Gad-el-Hak and Ho 9 also covered a 
wide range of k and observed the formation of additional 
vortices due to the motion of the trailing edge. They con- 
cluded that the complex flowfield of an oscillating airfoil 
can be explained primarily from the mutual induction be- 
tween the DSV and the "trailing edge shedding vortex". 

It also becomes apparent from a literature survey 
that the primary purpose of most previous experiments on 
dynamic stall was to determine the unsteady forces and 
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pitching moments on the airfoil 3,4 . Efforts to quantify the 
dynamic stall vortex, its evolution and trajectory, and its 
dependence on the reduced frequency k have been very 
limited. There have been several computational studies of 
the flowfield over a pitching airfoil (e. g. Shida et al. 10 , 
Visbal and Shang 11 ). These studies, however, require enor- 
mous computational power to obtain grid-independent solu- 
tion. Moreover, applicability and appropriateness of turbu- 
lence models required in the computation of the unsteady, 
separated flow under consideration are always questionable. 
Thus further experimental studies on the topic are clearly 
desirable. 

The mechanism for the generation and subsequent 
accumulation of vorticity leading to the formation of the 
dynamic stall vortex is discussed by Reynolds and Carr 12 , 
Ho 13 and Visbal 14 . However, validation of the proposed 
mechanisms is difficult due to a lack of detailed flowfield 
measurements. Also, most of the previous flowfield studies 
concentrated on the DSV while it is located over the airfoil 
surface, and the structure of the downstream wake remains 
relatively unknown. In only a few experiments was the 
vorticity distribution in the wake measured, but these exper- 
iments were limited in scope 15,16 . An understanding of the 
wake structure originating from the dynamic stall process is 
essential for analysis of more complex flows, e.g., in ma- 
chinery involving rows of blades 17 . These are the general 
issues which motivated the present experimental investiga- 
tion. 

The objectives of this investigation are twofold. 
First, to study the detailed flowfield for a specific k, and 
secondly, to study the effect of k on the dynamic stall pro- 
cess. The oscillation amplitude and the mean angle of attack 
are kept constant. Detailed flowfield measurements are 
conducted by the phase averaging technique. The distribu- 
tion of the spanwise component of vorticity is measured by 
phase averaging of the streamwise and the transverse com- 
ponents of the velocity. The measurements are performed 
for the actual spatial distribution over the cross-section of 
the flow, and thus, without invoking the Taylor hypothe- 
sis 18 . The unsteady vorticity data are analyzed to estimate 
the unsteady lift force on the airfoil, the concept of which 
is discussed in the text. Limited quantitative measurements 
are performed for other values of k over the range, O^k- 
^1.6. For all operating conditions, flow visualization is 
performed by the smoke-wire technique. 

2. Experimental procedure 

The experiments were carried out in a low-speed, 
76x51 cm wind tunnel as shown in figure la. The flow 
enters through a 16: 1 contraction section with five screens, 
passes through the test section and is then exhausted by an 
axial flow fan. The free-stream turbulence intensity in the 
test section is less than 0.1%. An NACA0012 airfoil with 
10.2 cm chord and 76 cm span is mounted horizontally at 


the center (mid-height) of the test section. The maximum 
blockage to the flow occuring at the maximum angle of 
attack of 25° is 8.5%. The airfoil is supported at the two 
ends by two 0.635 cm diameter rods each of which passes 
through a pair of cylindrical bearings housed in the tunnel 
walls. The pitching mechanism (figure lb) essentially con- 
sists of a crank and connecting-rod, together with a fly- 
wheel, to oscillate a lever arm. The lever arm (not visible 
in figure lb) oscillates an output shaft which in turn is 
connected to one of the airfoil support rods via a flexible 
coupling. The oscillation amplitude is adjusted by changing 
the crank radius and the oscillation frequency is controlled 
by varying the motor rpm. The entire pitching mechanism 
is installed in a steel frame which is structurally divorced 
from the tunnel and is secured to the vibration isolated bed 
plate of the test cell. 

For all data presented here, the airfoil is supported 
at the quarter chord location and the angle of attack is var- 
ied as a = 15° + 10° sin(2nft). In the following, the 
suffix "u" is used to indicate upstroke when a is increasing 
and H d" is used to indicate downstroke when a is decreas- 
ing. Experiments were conducted at two chord Reynolds 
numbers, R c = 22,000 and 44,000, and for k values of 0 
(fixed), .05, .1, .2, .3, .4, .6, .8, 1., 1.2 and 1.6. Experi- 
ments at the higher R c were conducted only for k^O.4. De- 
tailed measurements and analysis, however, are done for 
only one case, k = 0.2 (4 Hz) and R c = 44,000. 

The measurements are made using a single or a 
crossed hot-wire probe mounted on the traversing mecha- 
nism shown in figure la. The probe can be traversed in 
the streamwise (x) direction through a longitudinal slot in 
the test section floor. For a given x, the probe can be 
moved up and down in y through an automated computer 
control. The co-ordinates x and y are referenced to the 
airfoil trailing edge at 0° angle of attack, and unless 
otherwise stated all measurements reported are for the x-y 
plane at the mid-span location. The probe movement, 
data acquisition, and analysis were done by a MicroVAX 
3300 computer. 

The smoke wire technique is used for flow visu- 
alization. Essentially the same procedure was used as 
followed by Panda 19 . A 0.005 inch nichrome wire was 
placed vertically either upstream of the airfoil or just 
downstream of the trailing edge. A flash unit was used to 
illuminate the smoke streaks from the wire which were 
photographed using a 35mm camera. Control on the 
smoke-wire operation was not fine enough to produce 
photographs at a desired a during the oscillation. A small 
rod, marked every .25 inch by white dots, was placed 
downstream of the trailing edge. Photographs of the 
marker and the airfoil fixed at different a were taken to 
calibrate the marker. Usually, a large number of photo- 
graphs were taken at each flow condition. The position of 
the airfoil trailing edge relative to the marker determined 
the approximate angle of attack. The direction of motion 


was inferred from reviewing a large ensemble of photo- 
graphs. 

3. Experimental results 

Since the flow visualization photographs present 
the global view of the flowfield, they are discussed first. 
This is followed by the detailed measurements for the k 
= 0.2 case. The calculation of the unsteady lift from the 
vorticity flux is discussed next. Limited quantitative 
measurements on the effect of the reduced frequency k is 
discussed in the last section. 

3.1 Flow visualization 

3.1.1 Flowfield for k = 0.2. R. = 44.000 

Figure 2a shows a sequence of photographs at 
various phases of the oscillation cycle. The flow is from 
left to right. The marker with white dots is visible at the 
airfoil trailing edge. Frames A to F show phases when 
the angle of attack (a) is increasing (upstroke) and frames 
G to J show phases when a is decreasing (downstroke). 
As a increases, a clockwise vortex forms on the airfoil 
surface (frame D). This is the "dynamic stall vortex" as 
referred to by previous researchers. With further increase 
in a, the DSV moves towards the trailing edge. When it 
reaches the trailing edge, a counter-clockwise vortex 
starts to form near the trailing edge (frame F). This vor- 
tex becomes clearer in frame G and can be seen more 
clearly in figures presented later. The counter-clockwise 
vortex will be referred to as the trailing edge vortex 
(TEV) in the following. The TEV grows quickly beneath 
the DSV (frames F and G) and lifts the latter from the 
airfoil upper surface. The DSV and the TEV combine to 
form a structure whose cross section looks like a mush- 
room. This evolves and increases in size as it moves 
downstream (Frames H, I and J). In frame I, this vortical 
structure is about 2 X A chord lengths from the trailing 
edge. The transverse extent is already very large and 
measures about 3 chord lengths. 

The presence of the TEV can be seen in some 
flow visualization photographs and computational results 
of previous researchers. Gad-el-Hak and Ho 8 called it the 
trailing edge shedding vortex. Flow visualization for 
single pitch-up motion (pitch rate, a + = 0.1) by Walker 
and Robinson 20 showed such a vortex, which they ref- 
erred to as the "secondary vortex". Shida et al. 10 showed 
density contours, obtained computationally, at various 
phases for the flow around an airfoil oscillating between 
0° and 20° angles of attack at k = 0.1. The results 
showed an additional region of strong density gradient 
near the trailing edge at a = 19.3° when apparently the 
DSV was being shed. This can be interpreted as a conse- 
quence of the TEV. It will be shown later that the forma- 
tion of the TEV causes a large lift oscillation observed in 


the lower branch of the hystersis loop in the C f versus a 
curve. 

The ’mushroom-like’ structure observed in the 
wake has remained practically unnoticed in previous 
work. This enormous vortical structure, which occurs 
periodically at the oscillation frequency, may be quite 
significant in blade vortex interaction and aerodynamic 
noise generation from configurations involving rows of 
blades. After the passage of the DSV and the TEV, 
frames J and A indicate the passage of a few smaller 
vortices before another DSV begins to form on the upper 
surface. 

Figure 2b shows a sequence of smoke-wire pho- 
tographs for the same flow condition as in figure 2a but 
with the smoke-wire located near the trailing edge. 

Frames a, b and c represent phases when a is increasing 
and frames d and e represent phases when a is decreas- 
ing. Frames a,b,c,d,e in figure 2b approximately corre- 
spond to frames A,B,E,G and I in figure 2a, respectively. 
When the smoke-wire is kept upstream of the leading 
edge, the separated flow region on the upper surface and 
in the near wake remains dark due to the lack of convec- 
tion of the smoke-containing fluid. However, when the 
smoke-wire is placed downstream of the trailing edge, 
strong reverse flow at certain phases carries the smoke 
upstream, highlighting the reverse flow zones. The pho- 
tographs of figure 2b clearly show the extent of flow 
reversal at various phases during the oscillation cycle. 

Unlike a steady airfoil (see photographs later), 
strong reverse flows can be observed even at the mini- 
mum a = 5° for the oscillatory case (figure 2b). A com- 
putational study of the flowfield for an airfoil undergoing 
a single pitch-up motion by Visbal 14 also showed the 
presence of a strong reverse flow in a thin layer adjacent 
to the upper surface during the formation of the dynamic 
stall vortex. Once the separated flow reattaches and the 
DSV starts to form (frames B and b, figures 2a and 2b), 
no reverse flow is visible. Presumably, the flow is at- 
tached to the airfoil upper surface downstream of the 
DSV and thus the smoke is not entrained upstream. How- 
ever, it should be obvious that the DSV is accompanied 
by regions of reverse flow adjacent to the airfoil surface. 

3.1.2 Effect of the reduced frequency k 
Flowfield for k = 0, different a: Figure 3 shows visual- 
ization photographs of the fixed airfoil for three indicated 
values of a . The pictures in the right column are taken 
with the smoke-wire placed upstream of the leading edge 
and those on the left column are taken with the smoke- 
wire placed downstream of the trailing edge. The NACA- 
0012 airfoil undergoes complete stall at about a = 14°. 
The photographs for a = 15° and 25° show fully stalled 
flow. Due to the low operating Reynolds number, it is 
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apparent that a "laminar separation” has taken place even 
when the airfoil is at 5° angle of attack. It is interesting 
to note that with the smoke wire placed at the trailing 
edge, the separated region remains dark for a = 5°. This 
indicates weak reverse flow under laminar separation. Ex- 
periments on laminar separation also showed weak re- 
verse flows in such a situation 21 . Referring back to figure 
2b, note the difference with the oscillating airfoil case 
where at the same a there is a strong reverse flow. 

At larger angles of attack the airfoil wake is 
characterized by the Karman vortex street (alternatively 
referred to as the bluff-body shedding) which contains 
alternate clockwise and counter-clockwise vortices. The 
corresponding Strouhal number St # = f^Sina/U. was 
measured (f„ being the shedding frequency), to be 0.18, 
via spectral analysis of the u-signal. This is consistent 
with pervious observations 22 . 

Flowfield at different k: Figure 4 shows the flowfields for 
R c = 22,000 at different reduced frequencies, k. These 
photographs represent the approximate phases when the 
dynamic stall vortex is lifted up and the trailing edge 
vortex has clearly formed. At the lowest k ( = 0.05) the 
DSV and the TEV are not clearly visible and the flow- 
field appears to be in a quasi-steady state. As k is incr- 
eased above 0.1, the unsteady effects become prominent. 
The initiation of the "mushroom" structure is clear for the 
higher k cases. 

An interesting effect of varying k is the resulting 
phase delay in the shedding of the DSV. At k = 0.2, the 
shedding of the DSV takes place when the airfoil is near- 
ly at the end of the upstroke, at a = 24.5°. However, as 
k is increased, the shedding is delayed to the downstroke 
when a is decreasing. At k = 0.8, the DSV is shed 
when the airfoil is at the minimum angle of attack of 5°. 
A further increase in k brings the shedding of the DSV 
back to the upstroke part of the cycle. Similar observa- 
tions were made by McAlister et al. 23 , Kim and Seung 24 
and Leishman 4 in the limited k range of 0.05 to 0.25. 
With increasing k, a delay in the formation of the DSV 
was observed in the last experiment. In the experiment by 
Kim and Seung, a similar delay was observed in the 
formation of the separation region near the trailing edge. 
McAlister et al.’s data are similar to the present data and 
are further discussed in the following in connection with 
figure 5. 

From photographs similar to those in figure 2a, 
for k = 0.2, it is observed that the DSV starts to form at 
about a = 8° during the upstroke. In this case, the time 
elapsed until the airfoil reaches a = 25°, is enough for 
the DSV to move to the trailing edge and shed into the 
wake. However, when k is increased, the time needed by 
the DSV to form, grow and move to the trailing edge 
becomes larger compared to the time taken by the airfoil 


to complete the upstroke. This apparently causes the 
observed phase delay. At k = 0.8, the time needed by 
the DSV from formation (at the beginning of upstroke) to 
shedding equals the oscillation period and thus, the shed- 
ding occurs at the end of the downstroke. A further in- 
crease in k, to say 1.6, causes two DSVs: one that has 
started to grow at the leading edge and one left from the 
previous cycle near the trailing edge, to reside on the air- 
foil suction surface simultaneously. 

Figure 5 shows a plot of the angles of attack and 
the phase angle (o>t) at which the center of the DSV is 
approximately above the trailing edge before it sheds into 
the wake at various k. The values of a and a>t are related 
by a = 15° + 10°sin(a)t). Each data point is an estimate 
from several photographs similar to those shown in figure 
4. The phase data are cross checked with < U > con- 
tours, measured near the trailing edge as discussed later. 
The solid line in figure 5 represents data taken from 
McAlister et al. 23 for the same oscillation condition but at 
a much higher Reynolds number (2.5xl0 6 ). The latter 
data represent the instants when the suction peak near the 
leading edge collapses. The good agreement of these data 
with the present result, within the k-range covered, indi- 
cates that the phase delay phenomenon for comparable 
airfoils may be essentially independent of the Reynolds 
number. However, McAlister et al. also showed that 
leading edge modification such as by a boundary layer 
trip can significantly alter the phase delay characteristics. 

The delay in shedding of the DSV may have a signif- 
icant effect on the lift and pitching moments experienced 
by the airfoil at various k. It is known that the lift aug- 
mentation takes place with the formation of the DSV and 
the sudden drop of lift occurs when it leaves the airfoil 
surface. For low k, maximum lift as well as the sudden 
drop in lift occurs around the maximum angle of attack, 
a = 25°. At higher k, the shape of the hystersis loop in 
the C ( versus a curve is likely to change. For example, 
at k = 0.8, the maximum lift may occur around a = 5°. 
For still higher k, when two or more vortices exist on the 
upper surface simultaneously, complex hystersis loops 
may be expected in the C f versus a curves. 

Figure 6 shows pictures similar to those in figure 
4 but for different phases. The full grown DSVs (left 
column), and the "mushroom vortices" (right column), 
are shown for a few values of k. Again, neither the DSV 
nor the mushroom vortex are clear for k = 0.05. As k is 
increased, the DSV becomes stronger and the vortical 
motion in the core becomes intense. An expected trend 
with increasing k is the decrease in the spacing (wave- 
length) of the mushroom vortices in the wake. For k = 
0.2, only one mushroom vortex is seen in each frame 
covering a streamwise distance of about 3 chords, while 
for k = 1.6, spacing between two mushroom vortices 
decreases to 1.5 to 2 chords. 
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At k = 0.2, for which detailed flow visualiza- 
tion was shown earlier, the sequence of vortex formation 
is such that the DSV and the TEV together leave the 
trailing edge, and this is followed by the shedding of 
smaller vortices. As k is increased to 0.4 or higher, the 
scenario becomes somewhat different. A small clockwise 
vortex rolls up on the upper surface and is shed in the 
wake before the dynamic stall vortex; all other smaller 
vortices are suppressed. McAlister and Carr 25 observed 
this smaller vortex preceding the DSV in their flow 
visualization experiment and called it the "shear layer 
vortex". The present visualization pictures also appear to 
show that this vortex is accompanied by another counter- 
clockwise vortex prior to the shedding of the DSV. 

Note that at even higher k, there is a significant 
change in the "mushroom" structure. From figure 6 it can 
be seen that while the structure appears upright at k = 
0.6, it is tilted for k = 1.2, and appears upside down at 
k = 1.6. The vortex structures for k = 0.6 and 1.6 are 
sketched in figure 7. At the lower k, the DSV is shed 
before the TEV and both of them move upwards in posi- 
tive y. On the other hand, at k = 1.6, the TEV is shed 
earlier than the DSV and the resulting "inverted mush- 
room" remains in the same horizontal plane containing 
the pitch axis. A somewhat intermediate structure is seen 
for k = 1.2. 

A question may be asked if there exists an opti- 
mum range of k for which the vortices due to dynamic 
stall are the strongest and the largest. The visualization 
photographs shown in the preceding figures indicate that, 
within the distance covered in these photographs, the 
DSV and the TEV grow to the largest size for values of k 
in the range of 0.2 - 0.4 (see e.g., frame I, figure 2a). 

As k is increased further, the DSV as well as the "mush- 
room vortex" reduce in size. However, the vortical struc- 
tures become more complex at larger k. 

3.2 Flowfield measurements for k = 0.2, R c = 44,000 

3.2.1 Phase averaged vorticity measurements 

The axial and transverse velocity components, 
ensemble averaged over many oscillation cycles, are 
measured in the wake as well as on the suction side of the 
airfoil, avoiding the region covered by the pitching mo- 
tion. The periodic phase reference, once in every oscilla- 
tion cycle, is obtained by an optical pick-up mounted on 
the pitching mechanism (figure lb). At first, the signals 
from the crossed hot-wire probe and the optical pick-up 
are stored in a digital computer. The data are acquired at 
several x-stations. At each x-station the automated tra- 
verse is used to collect data at several transverse locations 
(y). These data are post-processed to obtain the phase- 
averaged distributions of < U > and < V > . Typically 
the averaging is done over 80 cycles or more. Finally, 
the spanwise vorticity component is calculated as, 


<^>(x 9 y,z) 


3<V> d<U> 
dx dy 


(1) 


The vorticity values are non-dimensionalized by UJc. 

Figure 8 shows examples of the phase averaged 

< U > profiles at three axial locations for five oscillation 
phases. For a given x, the instantaneous "wake profiles" 
vary dramatically with oscillation phase. At x/c = -0.3, 
the amplitudes are likely to be contaminated by hot-wire 
error due to reverse flow close to the airfoil surface (y/c 

< 0.15); the errors are discussed further in the follow- 
ing. The highest <U> variation at x/c = 0.6 is ob- 
served around the maximum a of 25° when the DSV and 
the TEV pass through the measurement station. The pro- 
file showing maximum variation of <U> shifts, due to 
the convection time taken by the mushroom structure, to 
a later phase at the downstream station of x/c = 1.5. Ve- 
locity profiles similar to those shown in figure 8, as well 
as the corresponding < V > -profiles, were obtained for 
50 phase points covering a complete cycle and 32 x-sta- 
tions. 

Equi-velocity contours of phase averaged < U > 
(left column) and < V > (right column) are shown in 
figure 9; the instantaneous values of a are indicated. The 
measurement domain is divided into two parts. The first 
region lies on the airfoil upper side: -0.75^x/c^0. 1 and 
O.l^y/c^l.O and is divided into 13 (y) x 12 (x) measur- 
ing stations; the second region covers the wake from near 
the trailing edge, .06^x/c^2.0 6 and -.75^y/c^l.0, and 
is divided into 20 (y) X 17 (x) stations. Since results 
from these two sets are patched together in each frame of 
figure 9, there is some mismatch at the interface. The 
mismatch is caused partly by measurement error and 
mostly by a slight phase drift between the two sets of 
measurements. As will be clear from the vorticity plots 
(shown shortly), the top row of figure 9 corresponds to 
an instant when large vortices are absent in the measure- 
ment domain. The middle row corresponds to a phase 
when the DSV is moving towards the trailing edge on the 
airfoil upper surface. The DSV is characterized by strong 
axial and transverse velocities around its core. The bot- 
tom row shows the velocity field when the "mushroom 
vortex" structure occupies the wake. These data clearly 
show the enormous changes in the unsteady flowfield that 
takes place periodically under the action of the vortical 
structures. The peak-to-peak variations in <U> and 

< V > can be observed to be comparable in these data. 

As indicated before, a crossed hot-wire probe was 
used to measure the velocity components, following simi- 
lar procedure as used by Westphal and Mehta 26 . There are 
some errors in these measurements, especially due to 
hot-wire rectification during instants of large flow angu- 
larity and flow reversal. Even without rectification the 
errors are significant when the instantaneous velocity 
vector lies outside ± 30° of the direction bisecting the 
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probe wire angle. There is no easy way to assess the 
errors accurately and the exact extent of the errors has 
remained unknown. The flow visualization photographs 
have shown that flow reversal occurs close to the airfoil 
upper surface when the DSV and the other vortices pass 
by. The measurement domain, however, excludes the 
regions near to the airfoil surface. Nevertheless, some 
flow reversal and large flow angularity are likely to have 
contaminated the data especially in the middle of the 
wake just downstream of the trailing edge. Contours of 
the flow angle (tan 1 <V>/<U>) were constructed from 
the measured <U> and <V> data (not shown). Down- 
stream of x/c = 0.25, the phase averaged flow angle was 
found to be confined in the range of -26.5° to 4-31.5°, 
for all measurement points at all phases. Inspection of the 
contours (which look very similar to the < V > -contours) 
also indicates that the data are free from rectification. On 
this basis, it is inferred that the errors due to large flow 
angularity are not significant downstream of one quarter 
chord from the trailing edge, as well as for y/c >0.15 
for the upstream measurement region above the airfoil. 
Good comparison of the measured vorticity field (dis- 
cussed shortly) with the flow visualization photographs 
supports this conclusion. 

3.2.2 Vorticity evolution 

The evolution of the spanwise component of 
vorticity at various phases is shown in figure 10. For k 
= 0.2, the distance through which the flow structures are 
convected downstream in one period of oscillation is 
much longer than that covered by the measurement do- 
main. Hence, each frame shows the wake structure for a 
small fraction of one complete cycle. The data acquisition 
rate was chosen such that the phase-averaged distributions 
are computed for 50 time instances per cycle. However, 
only 16 time instances, clearly capturing the passage of 
the vortices, are shown in figure 10. Note that the frames 
are not at equal time intervals. 

During the initial part of the cycle, as a is in- 
creased, the measurement region is relatively quiet and 
large vortices are absent. The wake region is marked by 
clockwise vorticity above the trailing edge, and counter- 
clockwise vorticity below it. As the angle of attack is in- 
creased above 8°, a small vortical region nucleates near 
the airfoil quarter chord and grows rapidly with accumu- 
lation of positive vorticity. Note that the equi-vorticity 
contours are drawn up to a maximum level of 12. Peak 
levels exceeding this value have been marked in some of 
the regions. The nucleation of positive vorticity over the 
airfoil is accompanied by an interesting change in the 
wake vorticity distribution. The negative vorticity in the 
lower part of the wake remains as is, but the positive 
vorticity in the upper part of the wake depletes dramati- 
cally (see, e.g., data for a = 20. 7u, figure 10). The 
positive vorticity, generated on the airfoil suction surface, 


which otherwise would have shed into the wake, starts to 
accumulate leading to the formation of the dynamic stall 
vortex. A similar process is described by Visbal 14 from 
his computational results. 

Increasing a above about 18° causes a large 
portion of the accumulated vorticity to move towards the 
trailing edge. This represents the dynamic stall vortex and 
is marked as VI in figure 10. The shedding of the DSV 
in the wake (around the phase marked 24. 9d), is accom- 
panied by the appearance of the counter-clockwise trailing 
edge vortex. Initially the TEV, marked V2, is small and 
lies underneath the DSV. However, as the two evolve, 
the TEV moves upward. The combination of the DSV 
and the TEV gives the appearance of a mushroom in the 
flow visualization photographs. During the downstroke, 
positive vorticity is shed continuously in the wake which 
rolls-up to form vortices V3 and V5. It is interesting to 
note that each of the clockwise vortices is followed up by 
a counter-clockwise vortex (V4, V6). Therefore, the 
wake contains a series of clockwise and counter-clockwise 
vortices after the passage of the DSV and the TEV. The 
succession of a counter-clockwise vortex following the 
passage of a clockwise vortex is, perhaps, expected as the 
net vorticity shed into the wake over a complete oscilla- 
tion cycle should be zero. The measured vortical struc- 
tures, their streamwise and transverse extents, agree quite 
well with the flow visualization photographs of figures 2a 
and 2b.‘ 

3.2.3 2-dimensionalitv of the flowfield 

In order to check the spanwise uniformity of the 
flowfield, phase averaged velocity measurements were 
conducted at various spanwise (z) stations for constant x 
and y. Figure 11 shows a set of <U> traces at different 
z/c, for x/c =1.0 and y/c = 0. More than one period of 
oscillation is shown in every trace. Passage of the dynam- 
ic stall vortex is marked by large positive peaks in the 
< U > traces. The velocity traces are quite similar and it 
is clear that the dynamic stall vortex arrives nearly at the 
same instant at all spanwise stations. Thus, the flowfield 
under study may be considered as reasonably 2-dimen- 
sional. 

3.3 Estimation of unsteady lift from the vorticity data 

An effort was made to estimate the time depen- 
dent component of the lift from the measured vorticity 
data. The variation of the unsteady component of lift over 
an oscillation period may be compared with the hystersis 
loop in C ( variation with a measured by others 27 . The 
underlying idea and the results are discussed in the fol- 
lowing. 

Let us consider the control volume ABCD 
around the airfoil as shown in figure 12. It is reasonable 
to assume that the flows at the boundaries DA, AB and 
BC are irrotational and all vortical fluid passes through 
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boundary CD. If the airfoil is held fixed at a given a, 
then the circulation around the path ABCD, T 0 , is a 
constant and the lift Lq is obtained via the Kutta-Joukow- 
skii theorem, L<, = pU So- For the steady flow, T 0 is 
equal in magnitude but opposite in sign to the total circu- 
lation shed in the starting vortex system. After the pas- 
sage of the starting vortex system, the circulation around 
the airfoil remains constant and the net vorticity flux 
passing through the boundary CD at all times should be 
zero 28 . In other words, at any time after the passage of 
the starting vortex system, the amount of positive and 
negative vorticity passing through the boundary CD 
should be equal and cancel each other so that J CD co z dy 
= 0 . 

Let us now consider that the airfoil oscillates 
periodically with a time period T p . The circulation, T(t), 
around ABCD is now a function of time, varying periodi- 
cally with period T p . At any instant the lift L(t) may be 
approximated by the Kutta-Joukowskii theorem, L(t) = 
pUJ^t) (see farther discussion at the end of this sec- 
tion). The vorticity flux through CD at any instant is 
now non-zero, and the net vorticity passing over a time 
At should be equal and opposite to the change in the 
circulation AT. Note that the vorticity flux through CD 
at any instant corresponds to the change in the circulation 
A T around ABCD at the same instant. Therefore, the 
vorticity flux through CD represents the variation in the 
unsteady circulation around ABCD without any time lag. 
If we assume that in time At the vorticity crossing the 
boundary CD is contained within the region CDD’C* 
(figure 12) then, 

Ar = r(r + Ar)-r(f) = j u^dy (2) 

Where, /w, dxdy represents circulation around CDD’C’. 
Invoking the Taylor hypothesis and using a suitable aver- 
age convection velocity, U^, one can write, 

Ar = -U^tj^dy (3) 


Where, the vorticity flux J CD o> z dy is a function of time. 
Use of a constant convection velocity in the Taylor hy- 
pothesis should represent a good approximation of the 
dynamics of the large scale vortical structures 18 , and 
therefore, the above equation should be a good estimate 
of A r. Letting 5(0 = J c d <*> z dy it follows in the limit 
that 


— = -u 5(f) 

dt ooov 


r(') - 


(4) 


Thus if the integration of the vorticity flux is started at 
time tQ, the unsteady component of circulation T(t) - IXt©) 
will be given by the integral in the above equation for 
any time t. Obviously, the integral over a complete peri- 
od T p should be zero as the magnitude of circulation 
should come back to the same value as that at the start of 
the period, i.e., -I- T p ) = T( to). T^) represents 

circulation at time ^ and remains unknown. However, the 
integration of the vorticity flux, carried through different 
intervals provides the unsteady component of circulation 
over a period. This in turn provides the unsteady compo- 
nent of the lift L^t) as a function of time (and therefore 
as a function of a) over a period. 

In the present case, the vorticity measurements 
are not so reliable near the trailing edge due to hot-wire 
errors (see section 3.2.1). Thus the integration of the 
vorticity flux is carried out at a downstream location x, 
over the boundary EF (figure 12), where measurements 
are relatively free from such errors. This, however, intro- 
duced a phase lag for the calculated variation of the circu- 
lation which can be estimated by applying the Taylor 
hypothesis once again. If fit = fix/U^ represents the 
time taken by the flow to reach EF from CD then, 

m = + r( fo ) (5) 


Where, 5,(t) is the vorticity flux through EF. The un- 
steady lift component at any instant of time is given as, 

L u (f) = ptur( 0 -r(g) = -p u.ur'Xm 


The unsteady component of the lift coefficient, nor- 
malized by 0.5pU„ 2 c, is plotted as a function of a in 
figure 13a. In the present calculation, = 0.6U„ is 
assumed based on the average convection velocity of the 
"mushroom vortex". The calculation is started when a = 
5° at the beginning of the upstroke. The axial station 
where the integration is performed is at x/c = 0.3, where 
the hot-wire measurement errors should no longer be a 
serious problem. It is worth mentioning here that varying 
the calculation station within the range 0.2^x/c^l.5 
produced an essentially similar C t variation as shown in 
figure 13a. 

The loop in the calculated C e vs. a in figure 13a 
does not close completely over one complete period. The 
amount by which the loop fails to close, reflects the 
cumulative error over one complete period (50 phase 
points) in determining the vorticity distribution at x/c = 
0.3. In spite of this error, the overall features of the 
calculated C ( variation are remarkably similar to the mea- 
surements of McAlister et al. 27 reproduced in figure 13b. 
Note that while figure 13a shows the unsteady component 


7 


only, the data of figure 13b includes the steady compo- 
nent. 

As a is increased, the dynamic stall vortex forms and 
increases in strength. The slope of C f versus a curve in- 
creases as the steady stall limit is exceeded. A compari- 
son with the flow visualization photographs of figure 2 
and vorticity evolution of figure 10 shows that the sudden 
drop in lift occurs when the DSV approaches the trailing 
edge and is lifted up from the airfoil surface. The re- 
bounding of the lift near the highest angle of attack dur- 
ing the downstroke is due to the shedding of the trailing 
edge vortex that effectively increases the circulation T(t) 
once more. The present data are characterized by undula- 
tions in the lower branch of the hystersis loop. The pas- 
sage of the combined DSV and the TEV is followed by a 
few other smaller positive and negative vortices causing 
the observed undulations. 

In spite of the overall similarity, clearly differ- 
ences exist between the data of figures 13a and 13b. 

Some of the differences may be attributed to different ex- 
perimental conditions between the present experiment (R c 
= 4.4xlO\ k = 0.2) and that of McAlister et al. (R c = 
5xl0 5 , k = 0.15). It should be emphasized here that large 
differences can occur in the hysteresis loop especially for 
different k, as clearly shown by McAlister et al 27 . and by 
the wealth of data obtained in earlier experiments at 
NASA Ames (see, e.g., Carr et al 3 .). Other factors, e.g., 
Reynolds number, aspect ratio of the airfoil and leading 
edge modification can also produce some changes in the 
hysteresis loop. The data in figure 13b were chosen from 
the cited reference such that the experimental conditions 
were the closest to those in the present case. However, 
note that both k and R c are quite different, and therefore, 
the differences between the two sets of data may not be 
surprising. There is also another important difference 
between the present calculation and the force measure- 
ment data. The force exerted by the fluid on an oscil- 
lating airfoil can be viewed as a result of two contribu- 
tions. The first one is due to the circulation around the 
airfoil and the second one is an inertia force associated 
with the mass of the fluid displaced by the airfoil (Wu 28 ). 
The present calculation scheme approximates the unsteady 
part of the former, but the latter, which can be shown to 
be small for low k (Visbal 14 ), is not accounted for. 

A conspicuous difference between the present 
result (figure 13a) and that in figure 13b is the undula- 
tions in the lower branch of the C ( curve. As mentioned 
before, the rebounding of the lift in the present data is 
followed by the undulations due to the alternate positive 
and negative vortices succeeding the TEV. Such undula- 
tions have been observed in other experiments (e.g., 
Leishman 4 ) and also in certain computational studies. 
Figure 13c shows a set of computational result obtained 
by a two-dimensional Navier Stokes code (L. N. Sankar, 


private communication). For details of the computational 
procedure the reader may refer to Wu, Huff and Sankar 30 . 
Even though these data are for a much different R c and a 
different k, and the use of turbulence models in unsteady 
flow calculations is subject to question, the overall fea- 
tures and the undulations in the lower branch are very 
similar to the present measurements. 

An advantage of the present procedure for determining 
the unsteady component of the lift is that it does not 
require pressure distribution or force measurements. The 
data is obtained solely from wake surveys. The pressure 
measurements may suffer from averaging effect due to 
inadequate distribution of the pressure taps as well as 
possible uncertainties due to inadequate dynamic response 
of the transducers. 

Apart from the hot-wire errors, that apparently lead to 
an unclosed hystersis loop, perhaps the main questionable 
aspect in the present calculation is the applicability of the 
Kutta-Joukowskii theorem. If the DSV is lifted up from 
the airfoil surface but still remains in the control volume 
ABCD in figure 12, the actual lift on the airfoil may 
drop; but, the circulation around the circuit ABCD would 
remain the same. Thus, the lift variation may not be truly 
represented by the circulation around ABCD. Note also 
that the latter is only approximated in the present calcula- 
tion. Another source of error, related to the above prob- 
lem, arises from the application of the Taylor hypothesis 
and the accurate estimation of the phase lag between the 
calculated force and the actual force experienced by the 
airfoil. 

3.4 Wake velocity spectra: 

The effect of pitching at various oscillation fre- 
quencies on the wake velocity spectra, measured one 
chord downstream and 0.6 chord above the pitch axis, is 
shown in figure 14. When the airfoil is not oscillating (k 
= 0) and it is at the maximum angle of attack of 25°, a 
strong peak at a Strouhal number St s = f 9 cSina/U„ = 

0.18 (where the shedding frequency, f„ = 26 Hz) is ob- 
served in the velocity spectrum. This peak corresponds to 
the bluff body shedding as discussed in connection with 
figure 3. When the airfoil is oscillated at low k, the u’- 
spectrum contains two peaks: one at the oscillation fre- 
quency and a second one at the bluff body shedding fre- 
quency corresponding to the maximum a. Flow visualiza- 
tion photographs for low k, similar to those in figure 4, 
do indicate the existence of the Karman vortex street and 
confirm the quasi-steady behavior of the flowfield. 

Since the bluff body shedding frequency is a 
function of the length scale cSina, it may be expected 
that over a period when a varies the vortex shedding will 
take place over a wide range of frequencies. However, 
the velocity spectra show a single peak at fcSina/U. = 
0.18, where a corresponds to the maximum angle of 
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attack. Since the shed vortices are the strongest and the 
airfoil dwells longer at the maximum a, the flowfield is 
dominated by vortices shed at this a. As the reduced fre- 
quency is increased the amplitude of the spectral peaks 
corresponding to the bluff body shedding decreases. 

Above k = 0.2 the bluff body shedding is completely 
suppressed and multiple harmonics of the pitching fre- 
quency appear in the velocity spectra. Again, this trend in 
the spectral data is in reasonable agreement with the flow 
visualization pictures discussed earlier. Since, a further 
increase in the oscillation frequency does not change the 
overall characteristics of the velocity spectra, spectra for 
k > 0.2 are not shown. 

3.5 The unsteady wake at fixed x for different k: 

Since the actual vorticity measurement, via the 
spatial data acquisition technique, involves extensive time 
and effort, such measurements were performed for only a 
single value of k (= 0-2). In order to further study the 
effect of k, phase averaged measurements were performed 
for a fixed x-station (x/c = 0.12) at different k. The x- 
station close to the trailing edge was chosen to clearly 
capture the phase delay in the shedding of the DSV with 
increasing k, even though there is some hot wire error at 
this station in the middle of the wake (section 3.2.1). 
Figure 15 shows the contours of <U> for six different 
k. The abscissa in these plots represents the phase posi- 
tion (t) normalized by the period of oscillation (Tp). 

Let us, at first, consider the data for k = 0.2 for 
which detailed measurements and flow visualizations have 
been discussed. The appearance of the dynamic stall vor- 
tex is characterized by a region of large < U > around 
x/Tp = 0.35 which correspond approximately to the 
angle of attack of 23°. This, in time, is followed by a 
thinner wake structure until the minimum a of 5° is 
reached at x/Tp = 0.9. This trend remains qualitatively 
the same for the other k values. However, as k is in- 
creased, the appearance of the DSV at the measurement 
station is clearly delayed. At k = 0.2, the appearance of 
the DSV occurs at close to the maximum a (a = 23°) 
during the upstroke, while at k = 0.6 this occurs close to 
the minimum a (a = 7°) in the downstroke. A further 
increase in k to 1.2 results in the DSV appearing back in 
the downstroke at a = 18°. These data confirm the 
phase lag in the shedding of the DSV observed in the 
flow visualization pictures. The contours of figure 15 also 
indicate that the vortex induced velocities are the stron- 
gest around k = 0.4 which falls in the range of k where 
the dynamic stall phenomenon was inferred to be the stro- 
ngest from the flow visualization photographs (section 
3.1.2). 

4. Conclusions 

The following are the primary observations and 
inferences made in the present study. 


(1) In addition to the well known dynamic stall 
vortex (DSV) forming on the suction surface, an intense 
vortex of opposite sense is observed to form near the 
trailing edge (TEV) just when the DSV is shed. The 
combined DSV and the TEV grow to a very large size 
and take the shape of a "mushroom" as they convect 
away from the airfoil. The transverse extent of this "mus- 
hroom" structure is the largest, more than two chords at 
two chord downstream of the trailing edge, for the inter- 
mediate k range of 0.2 to 0.4. It has been known that the 
shedding of the DSV from the airfoil surface causes a 
sudden drop in lift which is followed by a rebounding of 
the lift in the lower branch of the hystersis loop of C e vs. 
a curve. The present measurements and analysis show 
that this rebounding is caused by the passage of the TEV. 

(2) The flowfield is documented in unprecedent- 
ed detail using the phase averaging technique. The mea- 
sured evolution of the vorticity field over a complete 
oscillation cycle agrees well with the flow visualization 
photographs. 

(3) The unsteady circulation around the airfoil 
and therefore, the unsteady component of the lift, is 
estimated from the vorticity data by calculating the cumu- 
lative vorticity flux shed into the wake over one complete 
period of oscillation. The calculated unsteady lift is in 
reasonable agreement with the hystersis loop in the C e 
versus a curve measured by McAlister et al 27 . The tech- 
nique is a novel one, providing the unsteady aerodynamic 
forces strictly from wake measurements, and could be 
quite useful in future experiments of this kind. 

(4) At very low reduced frequency (k^O.l), the 
oscillating airfoil behaves in a quasi-steady manner. Wake 
velocity spectra and flow visualization show the domi- 
nance of the bluff body shedding that would have occ- 
urred if the airfoil were held at the maximum angle of 
attack. As the oscillation frequency is increased, the bluff 
body shedding is gradually suppressed and the dynamic 
stall vortex and the trailing edge vortex become the pri- 
mary features of the flowfield. 

(5) For k = 0.2 the dynamic stall vortex is shed 
nearly at the maximum angle of attack of 25°, but as the 
reduced frequency is increased, the shedding of the DSV 
is progressively delayed. At k = 0.8 it is shed nearly at 
the minimum a of 5°, and at k = 1.6 two DSVs reside 
simultaneously on the suction surface during some part of 
the oscillation cycle. 
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Fig. 2(a) 
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Fig. 6 Flow visualization photographs at indicated values of k; R* = 22,000. 
Left column: phase when the DSV approaches the trailing edge; a from top to 
bottom are 13°u, 22°u, 19°d, 6°d, 15°u. Right column: phase when the 
"mushroom vortex" passes through the near wake; a from top to bottom are 
21 °u, 21°d, 8°d, 7°u, 21°d. 
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Fig. 10 Iso-vorticity contours for k = 0.2, R* = 44,000. Instantaneous a, with "u" and "d" < 
and downstroke, are indicated. Contour levels are at an increment of 0.5. Dashed lines indicate 
clockwise) vorticity. 
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Fig. 12 Control volume for calculation of unsteady lift. Fig. 13 C t vs. a. (a) Present data for unsteady C e ; k = 0.2, R c = 44,000, 

M = .019. (b) Total C t measured by McAlister et al. 27 k = 0.2, R c = 4.8x 
10 5 , M = .036.; (c) Total C, computed by Sankar 30 * =.151, R c =3.45xl0 6 , 
M = .283. 




18 


Fig. 14 Wake velocity spectra measured at x/c Fig. 15 Temporal distribution of <U>/U„ measured at x/c = 0.12 at various k; R*. 

1.0 and y/c = 0.6 for different k; R« = 44,000. — 22,000. Contour levels are at an increment of 0.12; a versus t shown at top. 

Ordinate is in arbitrary logarithmic scale. Successive 
traces are staggered by one and a half minor divisions. 
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